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Abstract 

We give a modern approach to the famous Cardano and Ferrari formulas in 
the algebraic equations with three and four degrees. Namely, we reconstruct these 
formulas from the point of view of superposition principle in quantum computation 
based on three and four level systems which are being developed by the author. 

We also present a problem on some relation between Galois theory and Qudit 
theory. 
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1 Introduction 



The method to obtain algebraic solutions in (algebraic) equations has a very long history, 
see pQ . To solve a quadratic equation is very easy. To solve a cubic equation is not so easy 
and has been given by Cardano. How to solve a quartic equation is comparatively hard and 
Ferrari has given such solutions. After that many mathematicians including Lagrange, 
Gauss et al made how to solve equations refined. However, (algebraic) equations with 
degrees of more than four have in general no algebraic solutions, which was proved by two 
young mathematicians Abel and Galois. 

By the way, we are studying a quantum computation based on multi-level systems of 
(laser-cooled) atoms. The three and four level systems are particularly interesting. In 
the theories the generalized Pauli and Walsh-Hadamard matrices (gates) play a central 
role. We don't mention the topics, see for example [2] |B] as an introduction. 

The aim of this note is to show that some fundamental techniques of the generalized 
Pauli and Walsh-Hadamard matrices are applicable to construct the Caldano and Ferrari 
formulas. The essential point is the principle of superposition which is the heart of 
Quantum Computation. Namely, we reconstruct the Caldano and Ferrari formulas from 
the view point of this principle. 

Our method is so simple that it must have been known in some context. However, we 
cannot check all related references because the history is too long! We could not find such 
references within our effort. Anyway, by this note many students will be able to master 
the Caldano and Ferrari formulas without much effort. 

2 Mathematical Preliminaries 

We summarize the properties of Pauli matrices and Walsh-Hadamard matrix, and next 
state the corresponding ones of generalized Pauli matrices and Walsh-Hadamard matrix 
within our necessity. 
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Let {<7i, <7 2 , cr 3 } be Pauli matrices : 



0"! 



1 

1 



-i 

1 



0"3 



1 
-1 



By (JTJ) 02 = 'io"iO"3, so that the essential elements of Pauli matrices are {<7i,o- 3 } and they 
satisfy 



<7i = a, 



a x = (7i, 03 = (7 3 ; <t 3 (Ti = -cricr 3 = e cxi<7 3 . 



(2) 



The Walsh-Hadamard matrix is defined by 



V2 



1 -1 



G 0(2) c 17(2). 



(3) 



This matrix (or transformation) is unitary and it plays a very important role in Quantum 
Computation. Moreover it is easy to realize it in Quantum Optics as shown in jS]. Let us 
list some important properties of W : 



W 2 = 1 2 , w ] = w = w~\ 

o x = Wa 3 W~\ 



(4) 
(5) 



The check is very easy. 

Next let us generalize Pauli matrices to higher dimensional cases. Let {Ex, E 3 } be the 
following matrices in M(n, C) 



Si 





1 
1 

1 



1 



V 



1 



1 



a 



a 



V 



a 



n-l 



(6) 



where a is a primitive root of unity a n = 1 ( o = e ™ ) . We note that 



a = a n -\ l + a + '-' + a"- 1 = 0. 



The two matrices {£1, S 3 } are generalizations of Pauli matrices {a±, <7 3 }, but they 
hermitian. Here we list some of their important properties : 

En v-" 1 1 ■ vt y^n— 1 v-tf 



v-m— 1 
J 3 — ^3 



S3S1 — a S1E3 . 



For n = 3 and n = 4 Ei and its powers are given respectively as 

/ n 1 \ / n 1 \ 



Si = 



1 

1 

1 



1 
1 

1 



and 



1 

1 

1 
1 



1 

1 

1 

1 



1 
1 



If we define a Vandermonde matrix W based on a as 

/ 



W = 



1 1 

1 a^ 1 



l 



<7 



1 

(j2(n-l) 



1 a 11 - 2 a 2 ( n - 2 ) 



1 o- n " 2 cr 2 ^ 2 ) 
1 a"- 1 a 2 ^" 1 ) 



1 

^(n-l) 2 
(„-l)(„-2) 



T 2(n-1) 
a n-l 

1 

T 2(n-1) 



\ 



.(n-l)(n-2) 



r("-l) 2 



1 




then it is not difficult to see 



W*W = WW^ = l n , Ei = WE 3 W^ = WE 3 W-\ 



(12) 



Since W corresponds to the Walsh-Hadamard matrix (f3*]l. so it may be possible to call 

W the generalized Walsh-Hadamard matrix. If we write W* = (w a b), then 

1 1 (lixi \ 

w a b = —^& ab = -^exp ( ab ) for < a, b < n — 1. 

y/n \ n J 

This is just the coefficient matrix of Discrete Fourier Transform (DFT) if n = 2 k for some 

k 6 N, see 0. 

For n = 3 and n = 4 W is given respectively as 

/ i i 1 \ ( \ 1 i \ 



v/3 



1 1 1 

1 a 2 a 
1 a a 2 



1 

71 



1 



and 



W 



1111 

1 a 3 a 2 
I a 2 I a 
1 a a 2 a 3 



a 

2 



i -l-i\/3 -l+n/3 
2 2 

i -l+-t\/3 -l-ty/3 
2 2 



/l 1 1 1^ 



/ 



1 -I 
1 -1 
1 I 



-1 i 
1 -1 
-1 -z 



(13) 



(14) 



We note that the generalized Pauli and Walsh-Hadamard matrices in three and four 
level systems can be constructed (by using Rabi oscillations of several types) in a quantum 
optical manner, see [3] and 0. 

3 Cardano Formula 

We consider an algebraic equation with three degrees 



x 3 + ax 2 + bx + c = 0. 



(15) 



Then by replacing x — > x — (a/3) we have 



x + px + q = 



with p = b — y an d q = c — y + %r. However, we usually set 



x 3 + 3px + q = 



(16) 



with 



p — , q = c 

y 3 9 ' H 3 



ab 2a 3 
27 " 



For this equation we have the famous solutions called Cardano formula : Three solutions 
are given by 

(17) 



x = u + v , (tu + a 2 v , (T 2 u + av , 



where uq and v o are primitive solutions of binomial equations 

o —q + V*? 2 + 4p 3 o —q — Vq 2 + 4p 3 
u = , v = 



2wi 

and a = e s , 



The solutions (|17|) look like (classical) superpositions of Uq and Vq. From here we 
reconstruct the solutions by using generalized Pauli and Walsh-Hadamard matrices (JBJ), 
fll3(l in three level systems, which makes the principle of superposition more clearly. 

Let xi, x 2 , x 3 be three solutions of (fTB^) and set 



A 



( 



V 



x x 



x 2 



\ 



X3 J 



x 3 + 3pXj + q = 0. 



Remind that x\ + x 2 + £3 = from IjlfiJ). so trA = 0. Then we have 



A 3 + 3pA + qE = 0, 



(19) 



where is a unit matrix. Now we consider superpositions of solutions Xj (j = 1,2,3). 
For that we operate W in ()13j) on (fT9"j) as the adjoint action (W< W = E) 



(W ] AW) 3 + 3p{W^AW) + qE = 0. 



(20) 
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Let us calculate AW : 

/ 



W ] AW = - 



X\ + x 2 + x 3 x 1 + o 2 x 2 + crx 3 Xi + ax 2 + o 2 x 3 » 



X\ + OX 2 + (T 2 X 3 X\ + X 2 + £3 Xi + <J 1 X2 + CTX3 

y x\ + <7 2 x 2 + crx 3 %\ + crx 2 + o" 2 x 3 x\ + x 2 + x 3 y 



% /3 a \ / 

a p 
^ P a j 



a: 



1 

1 

1 



f/3 



1 

1 

1 



By the way, 



(0, [3, a) = (xi, x 2 , x 3 ) 



1 1 1 

1 a 2 a 
1 a a 2 



Conversely, 



Oi, x 2 , x 3 ) = (0, P, a) < 



f / 

1 

3 

I V 



1 


1 


1 


1 


2 

a 


a 


1 


a 


a 2 



'11^ 



} = (0, P, a) 



J J 



1 a a 2 
1 a 2 a 



Therefore we have 

x\ = a + P, x 2 = cr 2 a + aft, x 3 = aa + a 2 p. 

Compare this with (JT7J). Let us solve (J5D|). 

= (aSi + pTslf + 3p(aEx + PT 2 ) + qE 
= (a 3 + p 3 + q)E + 3a(aP + p)£ a + 3p(ap + p)S 2 

Therefore 

a 3 + p 3 = - q 
a P — ~P 



Since P = —p/a, we have 



a 



3 V 



ex 



3 -"9 



a 6 + ga 3 - p 3 = 0. 



If we set t = a , then t + qt — p 6 = 0, the solutions are just 



t = -q + Vq 2 + Ap A -q - ^q 2 + Ap s 



Since a and (3 are symmetric in the equation (|24j) . a and /3 are primitive solutions of 
equations 

3 -q + yjq 2 + Ap 6 3 -g - v / ?+4p 3 
a = g ' ^ = 2 ' ^ ^ 

This is our derivation of Cardano formula. 



4 Ferrari Formula 

In this section we consider an algebraic equation with four degrees. By the same reasoning 
in the preceding section we have only to treat 

x 4 + px 2 + qx + r = (27) 

For this equation we have the famous solutions called Ferrari formula : Let Ao be a 
solution of the (decomposition) equation with three degrees 

A 3 - p\ 2 - 4rA + (Apr - q 2 ) = 0. (28) 

Then four solutions of (|27|) are given by ones of following two equations with two degrees 

* 2± v^{- ^} + ^°- (») 

This formula is not similar to that of Cardano. From this formula we cannot see the 
principle of superposition (of solutions), so we are dissatisfied at this one. From here we 
reconstruct the solutions by using generalized Pauli and Walsh-Hadamard matrices 0, 
()14|) in four level systems to make the principle of superposition clear. Our method may 
be fresh to not only students (non-experts) but also mathematicians. 
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Let x%, x 2 , x 3 , x 4 be four solutions of (}2"Tj) and set 



.4 



\ 



x 2 



x 3 



V 



Xj + px 2 + qxj + r = 0. 



(30) 



x 4 J 



Remind that x\ + x 2 + x 3 + x 4 = from (J2IJ), so XxA = 0. Then we have 

A* +pA 2 + qA + rE = 0, 



(31) 



where is a unit matrix. Similarly in the preceding section we consider superpositions 
of solutions Xj (j = 1,2,3,4). We operate W in ffT^j) on (}3~Tj) as the adjoint action 
(WW = E) 

[W ] AWY + p[W ] AW) 2 + qW ] AW + rE = 0. (32) 
The calculation of AW is as follows : 



W ] AW 
I 



\ 



a 7 
7 a /3 
^ 7 « j 
aEi + 7^ + /3£?, 



a 



1 

1 

1 
1 



+ 7 



1 

1 

1 

1 



1 

1 
1 

1 



(33) 



where 



I 1 

= -(xi + x 2 + x 3 + x 4 ), a = ~(xi + ax 2 + a 2 x 3 + a 2 x 4 ), 

7 = -(xi + a 2 x 2 + x 3 + a 2 x A ), (3 = ^(xi + a 3 x 2 + o 2 x 3 + cxx 4 ) 



(34) 



and a = e 4 = %, Then 



(0, P, 7, a) = (xi, x 2 , x 3 , x A ) 



'l 1 1 ^ 

1 a 3 a 2 a 

1 a 2 1 a 2 

1 a a 2 a 3 



Conversely, 

(xi, x 2 , x 3 , ar 4 ) 



= (0, /3, 7, «) < 



I V 



'l 1 1 1 M 



1 a 3 a 2 a 
1 a 2 1 a 2 
1 a a 2 a 3 



-i 



/J 



(0, /?, 7, a) 



V 



1 


1 


1 


1 


1 


a 


a 2 


a 3 


1 


a 2 


1 


a 2 


1 


a 3 


a 2 


a 



J 



(35) 



so we have 

27 = a + 7 + P, x 2 = cr 3 a + a 2 ^ + aft, 
£3 = a 2 a + 7 + a 2 p, £4 = era + cr 2 7 + er 3 /3. 

This is just superpositions of a, 7 and p. Next let us determine them. From (}3*3*j) 

(W ] AW) 2 = ( 7 2 + 2a/3)£ + 2 1 P^ l + (a 2 + p 2 )Y>\ + 2a 7 £ 3 

and 

(WUW) 4 = {{W ] AW) 2 } 2 ={(a 2 + P 2 ) 2 + ( 7 2 + 2aP) 2 + 8«/? 7 2 } E + 

4 {(« 2 + /3 2 )a + ( 7 2 + 2aP)p) 7 S! + 
2(a 2 + /3 2 )(3 7 2 + 2a/3)£ 2 + 
4 {(« 2 + /3 2 )/? + ( 7 2 + 2a/3)a} 7 £ 3 . 

Then the equation (J32*|l gives 

(a 2 + P 2 ) 2 + ( 7 2 + 2ap) 2 + 8a/? 7 2 + p( 7 2 + 2ap) + r = 
4 {(a 2 + p 2 )a + ( 7 2 + 2ap)p} 7 + 2p 7 /3 + qa = 
2(a 2 + /3 2 )(3 7 2 + 2a/3) + p(a 2 + p 2 ) + g 7 = 
^ 4 {(a 2 + p 2 )p + ( 7 2 + 2a/3)a} 7 + 2pa 7 + g/3 = 
Now let us transform the equations into a more convenient form (this point is important) 

(a) (a 2 + P 2 ) 2 + ( 7 2 + 2ap) 2 + 8a/? 7 2 + p( 7 2 + 2a/3) + r = 

(6) 2 7( 5(27 2 + 6ap + p) + 4a 3 7 + qa = 

(c) (a 2 + /3 2 ) (6 7 2 + Aap + p) + g 7 = 

(d) 2a 7 (2 7 2 + 6ap + p) + 4/3 3 7 + qp = 

10 



(36) 



(37) 



From (b) and (d) 



4a 3 7 + qa 



so 



2 7 P 

4a 4 7 + qa 2 = 4/3 4 7 + q{3 2 



2 7 2 + Qaf3 + p 



4/? 3 7 + q/3 



(a 



2«7 

-P 2 ) {A(a 2 + p 2 ) 1 + q} = 0, 



so we have 



a 2 + p 2 



Q_ 
47' 



(3* 



Substituting this into (c) 



P 



P 



67^ + 4a/3 + p-4 7 i = 0^7^ + 2a/3=-^ or 2a/3 = - 7 2 - (39) 



Substituting these into (a) 



PA 2 , / 2 . P\ 2 P 2 . n 



(40) 



Rearranging this equation 

64 7 6 + 32p7 4 - 4(4r - p 2 ) 7 2 - q 2 = 0. 
If we set T = 7 2 , then wa obtain the algebraic equation with three degrees 

64r 3 + 32pr 2 - 4(4r - p 2 )T - q 2 = 0. 

Since this equation can be solved (in the preceding section), we obtain 7 = yT for the 
largest real solution T. Next we look for a and p. From ()38j) and ()39jl 

« 2 + /9 2 = -^ 
a/3 = -|(r + f) 

By the same method in (|2*4^) a and P are primitive solutions of equations 

1 ' 1 ' 1 ' 1 ^ 1 (42) 

(43) 



(41) 




~ 2 I 4 7 ^ ^4 7 

Namely, we determined a, (3 and 7 completely. 

This is our derivation (a variation) of Ferrari formula 
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5 A (Real) Problem 



In this section let us consider an algebraic equation with five degrees. 



x 5 + px 3 + qx 2 + rx + s = 0. 



(44) 



We can apply the same method in the case of three and four degrees to this equation, 
however we cannot solve this one in an algebraic manner as in the preceding sections 
(Abel and Galois). We know the Galois theory to decide whether an algebraic equation 
has algebraic solutions or not. This is a classical theory. 

Then we have some natural questions : Is there some relation between the Galois 
theory and the principle of superposition ? What is a role of Galois theory in general 
Qudit theory ? Moreover, is it possible to consider a quantum Galois theory in general 
Qudit theory ? 

These are interesting problems in Quantum Computation. 
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Akira Asada pointed out to me that my method to derive the variation of Ferrari 
formula was similar to that of Euler which I didn't know when writing this note. In fact 
his point was correct, so let us copy the Euler formula in [iOj 1 for the completeness (this 
formula is not necessarily popular even in Mathematics). 

For the quartic equation (|27j) we consider the cubic equation 



Appendix Euler Formula 




(45) 



I could not find a textbook written in English which comments on the formula 
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We note that this equation is just our (reduction) equation (}4TJj) . Since this equation has 
three solutions (ti, t 2 , t 3 ), we choose signs of square roots y/ti, v^3 to satisfy 

VhVhVh = -~ . (46) 

o 

Then the solutions are given by 

xi = Vh. + Vh + Vh, x 2 = ~ Vh - Vh, 

x s = -Vh + Vh~ Vh, x 4 = -Vh- Vh + Vh- (47) 

To derive the solutions one use various relations between solutions (roots) and coefficients, 
which is quite complicated. 

Our method is just a modern derivation of the Cardano and Euler formulas which 
doesn't use relations between solutions (roots) and coefficients in an apparent manner. 
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